Astronomy & Astrophysics manuscript no. aa 


June 10, 2008 


(DOI: will be inserted by hand later) 





XMM- Newton studies of a massive cluster of galaxies: 

RXCJ2228.6+2036 

S.M. Jia 1,2 , H. Bohringer 1 , E. Pointecouteau 3 , Y. Chen 2 , and Y.Y. Zhang 4 

qq ■ 1 Max-Planck-Institut fur extraterrestrische Physik, 85748 Garching, Germany 

| e-mail: jiasm@ihep.ac.cn 

2 Key Laboratory of Particle Astrophysics, Institute of High Energy Physics, Chinese Academy of Sciences, 
(N 1 Beijing 100049, P.R. China 

3 CESR-CNRS, 9 Av. du Colonel Roche, 31028 Toulouse, France 

4 Argelander-Institut fur Astronomie, Rheinische Friedrich-Wilhelms-Universitat Bonn, Auf dem Huegel 71, D- 
53121 Bonn, Germany 



Oh 



Received / accepted 

Abstract. We present the X-ray properties of a massive cluster of galaxies (RXCJ2228. 6+2036 at z — 0.421) 
using XMM-Newton data. The X-ray mass modeling is based on the temperature and density distributions of 
the intracluster medium derived using a deprojection method. We found that RXCJ2228. 6+2036 is a hot cluster 
(Tsoo = 8.92]tj'32 keV) showing a cooling flow rate of 12.0^^ o Meyr -1 based on spectral fitting within the 
cooling flow radius (r coo ; = 147 ± 10 kpc). The total cluster mass is M 5 oo = (1.19 + 0.35) x 10 15 M and the mean 
gas mass fraction is f gas — 0.165 + 0.045 at rsoo = 1.61 + 0.16 Mpc. We discussed the PSF-correction effect on the 
spectral analysis and found that for the annular width we chose the PSF-corrected temperatures are consistent 
with those without PSF-correction. We observed a remarkable agreement between X-ray and SZ results, which is 
of prime importance for the future SZ survey. RXCJ2228. 6+2036 obeys the empirical scaling relations found in 
general massive galaxy clusters (e.g. S-T, M-T, L-T and M—Y) after accounting for self-similar evolution. 
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^5 1. Introduction that RXCJ2228. 6+2036 is a hot and massive cluster with 
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T = 10.4+ 1.8 keV, M tot (r < R v = r 17S = 2.9 Mpc) 
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^ The gravitational growth of fluctuations in the matter (1.8 ± 0.4) x 10 15 M Q , and with a gas mass fraction of 

density distribution can be traced by the evolution of the f gas = 0.22 ± 0.06 (Pointecouteau et al. 2002). Recently, 

galaxy cluster mass function (e.g. Schuecker et al. 2003). LaRoque et al. (2006) performed the Chandra X-ray ver- 

The hot and distant clusters are at the upper end of the sus OVRO/BIMA interferometric SZ effect measurements 

5-j ; mass distribution, thus they can be used to probe the cos- for the same cluster, giving T = 8.43±g;^ keV, f gas = 

mic evolution of large-scale structure and are therefore 0.138 + 0.009 from the X-ray data, and f gas = 0.1881°°^ 

fundamental probes for cosmology. But to date still very from the sz data at r250Q RXCJ2228.6+2036, as one of 

few hot and distant clusters are known. Therefore, it is the c i usters m the sample of X-ray luminous galaxy clus- 

important to study such clusters in detail, especially in terg with both x _ ray (Chandra) and SZ observations in 

X " ra y- Morandi et al. (2007), has a temperature of T = 6.86± ;|? 

RXCJ2228. 6+2036 is one of the distant (z = 0.421) keV and a total mass of M tot = (4.90 ± 4.35) x 10 14 M Q 

and X-ray luminous clusters of galaxies in the northern at r 5 oo = 1033 ± 464 kpc. However, the above results are 

sky. It is suspected to be massive and hot, and was well all based on the mass modeling under the assumption of 

recognized as an extended X-ray source in the ROSAT All- isothermality of the ICM. The XMM-Newton EPIC in- 

Sky Survey, included in both the NORAS galaxy cluster struments have both high spatial and spectral resolutions 

survey (Bohringer et al. 2000) and the ROSAT Brightest and a large field of view, and are therefore suitable for 

Cluster Sample (Ebeling et al. 2000). a spatially resolved spectral analysis. We make use of 

The first combined SZ versus X-ray analysis for XMM-Newton observations to carry on a detailed study 

RXCJ2228.5+2036 is based on the SZ data from the of RXCJ2228.6+2036 based on the X-ray mass modeling 

Nobeyama Radio Observatory (NRO) 45 m radio tele- usm S a spatially resolved radial temperature distribution 

scope and the X-ray data from ROSAT/HRI. It shows and Perform a detailed X-ray versus SZ comparison. 
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The structure of this paper is as follows: Sect. 2 de- 
scribes the data, background subtraction method and 
spectral deprojection technique. Sect. 3 presents the spec- 
tral measurements using different models to derive the ra- 
dial temperature profile, cooling time and mass deposition 
rate. In Sect. 4 we show the radial electron density profile 
and X-ray mass modeling. In Sect. 5 we discuss the im- 
pact of the PSF correction on the spectral analysis, and 
compare RXCJ2228. 6+2036 to the SZ measurements and 
the empirical scaling relations for massive galaxy clusters. 
We draw our conclusion in Sect. 6. 

Throughout this paper, unless explicitly stated other- 
wise, we use the 0.5-10 keV energy band in our spectral 
analysis. The cosmological model used is H = 70 km 
s _1 Mpc _1 , TO = 0.3, and fl\ = 0.7, in which 1' corre- 
sponds to 332.7 kpc at the distance of RXCJ2228.6+2036. 

2. Observation and data reduction 

2.1. Data preparation 

RXCJ2228. 6+2036 has been observed for 26 ksec in 
November 2003 by XMM-Newton and its observation ID 
is 0147890101. For our purpose, we only use EPIC data 
(MOS1, MOS2 and pn). The observations are performed 
with a thin filter and in the extended full frame mode for 
pn and the full frame mode for MOS. Throughout this 
analysis, we only use the events with FLAG=0, and with 
PATTERN<4 for pn and PATTERN<12 for MOS. The 
reduction was performed in SAS 7.1.0. 

The light curve of the observation shows some flares 
(i) in the hard band (above 10 keV for MOS and above 12 
keV for pn), possibly caused by the particle background, 
and (ii) in the soft band (0.3-10 keV), possibly due to the 
episodes of 'soft proton flares' (De Luca & Molcndi 2004). 
Therefore both the hard and the soft bands are used to 
select the good time intervals (GTI) as described in Zhang 
et al. (2006). The GTI screening procedure gives us 22 ks 
MOS1 data, 22 ks MOS2 data, and 18 ks pn data. 

We applied the SAS task edetect_chain to detect the 
point sources (the radius of the point sources is 0.6' 
containing 93% flux from the point source), and 
excised all of them from the cluster region. Then, a SAS 
command evigweight was used to create the vignetting 
weighted column in the event list to account for the vi- 
gnetting correction for the effective area. 

Due to read-out time delay, the pn data require a cor- 
rection for the Out-of-time (OOT) events. We created the 
simulated OOT event file and used it to account for this 
(see Striider et al. 2001 ) in our analysis. 

2.2. Background subtraction 

The XMM-Newton background normally consists of the 
following two components, i) Particle background: high 
energy particles such as cosmic-rays pass through the 
satellite and deposit a fraction of their energy on the de- 
tector. It dominates at high energies and shows no vi- 



gnetting, ii) Cosmic X-ray Background (CXB): the CXB 
varies across the sky (Snowden et al. 1997). It is more 
important at low energies and shows a vignetting effect. 

We choose the blank sky accumulations in the Chandra 
Deep Field South (CDFS) as the background (Zhang et al. 
2007), which was also observed with the thin filter. We ap- 
plied the same reduction procedure to the CDFS as to the 
cluster in the same detector coordinates, and the effective 
exposure time we obtained for CDFS is 54 ks for pn, 61 
ks for MOS1 and 61 ks for MOS2. RXCJ2228. 6+2036 is a 
distant cluster (z = 0.421), and we estimated that the 
signal-to-noise ratio of the region 6' < R < 6.5' is 
about 20%, so we can approximately assume that 
the emission of the cluster only covers the inner 
part of the field of view (R < 6'). The outer region 
(6.5' < R < 8') can thus be used to monitor the residual 
background. Here we applied a double-background sub- 
traction method to correct for these two kinds of back- 
ground components as used in Arnaud et al. (2002). First 
we estimate the ratio of the particle background, a, be- 
tween RXCJ2228. 6+2036 and CDFS from the total count 
rate in the high energy band (10-12 keV for MOS and 12- 
14 keV for pn), as described in Pointecouteau et al. (2004). 
So and Bo are the background spectra of the cluster and 
CDFS in the region of 6.5' < r < 8' with an area of A , 
and Si and Bi for spectra in the ith ring of the cluster 
and CDFS with an area of Ai. Then the cluster spectrum, 
after the double-background subtraction, S(i) is (e.g. Jia 
et al. 2006, Zhang et al. 2006): 

S(i) = S i -aB i - ^(S -aBo) (1) 
2.3. Spectral deprojection 

The combined image of MOS1 and MOS2 in the en- 
ergy band 0.5-10 keV is shown in Fig. 1. It is corrected 
for vignetting and smoothed with a maximum Gaussian 
smoothing size of a ~ 5 pixels. As shown in this figure, 
the X-ray emission of RXCJ2228. 6+2036 appears to be 
extended and almost symmetric except for some bright 
point sources (which were subtracted before the spectrum 
extraction, see Sect. 2.1). We extract the spectra from 
annular regions centered on the emission peak, and the 
width of each ring is determined according to the crite- 
rion described in Zhang et al. (2007): ^2000 net counts 
per bin in 2-7 keV to get a temperature with ^15% un- 
certainty. Considering the PSF (Point Spread Function) 
effect of XMM-Newton EPIC, whose Full Width at Half 
Maximum (FWHM) is 5" for MOS and 6" for pn, the min- 
imum width of each ring was set at 0.5'. We thus obtained 
5 annuli to extract spectra out to 6'. 

The deprojected spectra are calculated by subtracting 
all the contributions from the outer regions. Within each 
radial range, we assume the same spectrum per unit vol- 
ume. The deprojected spectrum of the ith shell is then cal- 
culated by subtracting the contributions from the i+lth 
shell to the outmost one (e.g. Matsushita et al. 2002 and 
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Fig. 1. The combined image of MOS1 and MOS2 of 
RXCJ2228. 6+2036 (12' x 12') corrected for vignetting and 
smoothed with a maximum Gaussian smoothing size of a = 5 
pixels. 



to 0.421 and the absorption Nh to the Galactic value 
4.68xl0 20 cm~ 2 (Dickey & Lockman 1990). The fitting 
results are listed in Table 1 and the central spectra fitted 
by this model are shown as Fig. 3(a). 

The deprojected temperature profile shows a drop in 
the core and a decrease in the outer regions (see the dia- 
monds in the upper panel of Fig. 2), which can be fitted 
by the following formula (Xue et al. 2004): 

n m , A i (lnr-lnr ) 2 

T(r) = T + —j— exp{ . 4) 

r/r to 

The best fit parameters are: A = 4.880 ± 0.001 keV, 
r = 2.494 ± 0.003', lu = 2.232 ± 0.004, T = 3.084 ± 0.001 
keV, xled = 0-31, and the best-fit profile is shown as the 
solid line in the upper panel of Fig. 2. From the tem- 
perature distribution, we can estimate the normalization- 
weighted temperature within 6', t mean = 8.57j^5 6 keV, 
which is consistent with the results of Pointecouteau et 
al. (2002) and LaRoque et al. (2006) within the error bars 
and a little higher than that of Morandi et al. (2007). The 
diamonds in the bottom panel of Fig. 2 show the depro- 
jected abundance distribution of RXCJ2228. 6+2036. 



Nulsen & Bohringer 1995). The detailed calculation pro- 
cedures are described as in Jia et al. (2004, 2006): 



D(i) 



N 



S(i)- C v (k,i)-D(k) 
fc-i+i 



/C v (i,i) 



(2) 



here D(i) is the deprojected spectrum of the ith shell, S(i) 
is the double-background subtracted spectrum of the ith 
shell and C v (k, i) is the fraction of the volume of the fcth 
shell projected to the ith ring. 

The on-axis rmf (response matrix file) and arf (auxil- 
iary responds file) are generated by the SAS task rmfgen 
and arfgen and are used to recover the correct spectral 
shape and normalization of the cluster emission compo- 
nents. 

3. Spectral analysis 

3.1. Radial deprojected temperature profile 

The spectral analysis is carried out in XSPEC version 
11.3.2 (Arnaud 1996). To study the temperature distri- 
bution of RXCJ2228. 6+2036, we perform a joint fit to the 
spectra of pn and MOS with an absorbed Mekal model: 



Modelx = Wabs(N H ) x Mekal(T, z, A, norm) 



(3) 



in which Wabs is a photoelectric absorption model 
(Morrison & McCammon 1983) and Mekal is a single- 
temperature plasma emission model (Mewe et al.1985, 
1986; Kaastra 1992; Liedahl et al. 1995). The tempera- 
ture T, metallicity A and normalization (emission mea- 
sure) norm are free parameters. We fixed the redshift z 




2 3 4 

Radius (arcmin) 




Fig. 2. Upper panel: Radial temperature profile of RXCJ 
2228.6+2036. Diamonds: the deprojected temperature and the 
solid line is the best-fit profile. Stars: for the PSF-corrected 
temperature (see Sect. 5.1). We offset the stars 5" to the left so 
as to illustrate these two kinds of temperature clearly. Bottom 
panel: Radial deprojected abundance of RXCJ2228. 6+2036. 
The confidence level is 90%. 



3.2. Mass deposition rate 

The temperature drop in the central part of RXCJ2228.6 
+2036 might indicate the existence of a cooling flow in 
the center. We thus estimate the parameters of the cooling 
flow as follows. 
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Table 1. The best-fit free parameters of RXCJ2228. 6+2036 by the single-temperature model: the temperature T; the abundance 
A and the normalized constant norm for the simultaneously fitting of pn and MOS. norm = 10~ 14 /(47t(Da x (1+z)) 2 ) j UenndV , 
where Da is the angular size distance to the source (cm) and n e is the electron density (cm -3 ). L bo i represents the bolometric 
luminosity (0.01-60 keV) in the units of 10 44 erg s _1 . The errors represent a confidence level of 90%. 



r(') 


T (keV) 


A (solar) 


noTO^lC^cm" 5 ) 


Xtod/dof 


Lf,oi(10 44 erg s" 1 ) 


0.0-0.5 
0.5-1.0 
1.0-1.75 
1.75-3.0 
3.0-6.0 


8.26iS;S 
11.65I?;? 4 

o Q1 +1.66 


31+ U iy 

u - ol -0.18 

0.20+°f 
0.30lJ;S 
12+ ' 25 
0.12(fixed) 


0.72 ± 0.03 

0.88 ±0.05 

1.08 ±0.06 

1.28 ±0.07 

7S+ 013 
u -'°-o.io 


1.12/251 
1.02/262 
1.00/196 
0.86/189 
0.51/59 


4 65 +,) ' 38 

a 70+O.58 
' o -0.76 
7 oq+0. 72 
' • oo -0.74 

7.89±; : ii 
04+1.20 



Table 2. The cooling time and the cooling flow rate deter- 
mined with the spatial method of the inner two regions of 
RXCJ2228.6±2036. The errors are at the 68% confidence level. 
1' = 332.7 kpc. 



ri 


r2 




MtMoyr- 1 ) 


0' 


0.5' 


1.18 ±0.07 x 10 iu 


200.2 ± 12.4 


0.5' 


1' 


3.58 ±0.27 x 10 10 


414.7 ±32.6 



The cooling time t coo i is the time scale during which 
the hot gas loses all its thermal energy, which is calculated 
as (e.g. Chen et al. 2007): 



tcool 



5 n e ± rii UbT 



(5) 



where A(T) is the cooling function of the gas, and n e , rii 
and riH are the number densities of the electrons, ions and 
hydrogen, respectively. Here for the nearly fully ionized 
plasma in clusters, n e = 1.2nu and n, = l.lnn- The 
determination of n e is explained later in Sect. 4.1. The 
cooling time t coo i of the inner two regions is given in Table 
2. The cooling radius designates the region inside which 
the hot gas loses all its thermal energy within a cluster life 
time scale, usually using the age of the universe (1.04 x 
10 10 yr at z = 0.421). The resulting cooling radius for n e {r) = n i 
RXCJ2228.6±2036 is r cool = 147 ± 10 kpc. 

We also fit the central spectra of pn and MOS by 
adding a standard cooling flow model to the isothermal 
Mckal component: 

Modeh = Wabs(N H ) x (Mekal(T, z, A,norm) + 

Zwabs(AN H ) x Mkcflow(M)). (6) 

Wabs and Mekal have been described in Sect. 3.1, Zwabs 
is an intrinsic photoelectric absorption model (Morrison & 
McCammon 1983), and Mkcflow is a cooling flow model 
(Fabian, 1988); ANh is the intrinsic absorption and M is 
the rate of gas cooling out of the flow. The fitting results 
show that the mass deposition rate in the central region 
r < 0.5' is M.OI^o M 0y r_1 ( see Table 3 )> so within r cool 
M = 12±f| M yr _1 . Fig. 3(b) presents the central spectra 
fitted by this model. 



4. Mass determination 

4.1. Electron density 

We divided the r < 6' region into 13 annuli centered on 
the emission peak, where the width of each annular region 
is determined to obtain at least 2000 total counts in each 
annulus region in 0.5-10 keV. After the vignetting correc- 
tion and the double-background subtraction, the surface 
brightness profile for RXCJ2228.6±2036, S(r), is derived, 
which can be fitted by a double-/? model (as Eq. 9) con- 
volved with the PSF matrices (Ghizzardi 2001) to correct 
for the PSF effect. Using the deprojection technique de- 
scribed in Sect. 2.3, we deproject the double-/? model (PSF 
corrected) and obtain the count rate of each correspond- 
ing shell, Ctr(i). Since the temperature and abundance 
profiles are known, giving T(i), A(i) and Ctr(i) in the ith 
shell, we can derive the corresponding norm(i) in XSPEC 
with Model\. The radial electron density n e of each region 
can be determined from Eq. (8), shown as stars in Fig. 4. 



norm(i) = 10" 14 /(4ttL> 2 ) • / n e n H dV. 



(7) 



We fit the derived electron density (the stars in Fig. 
4) with the double-/3 model (Chen et al. 2003): 



r 2 
l + ( — ) 



-n 2 



r 2 

l + ) 

r C 2 



,(8) 



the best-fit parameters are: noi = 0.0027 ± 0.0002 cm 3 , 
r cl = 2.5480 ± 0.0002 arcmin, fa = 1.4031 ± 0.0001, n 02 
= 0.0109 ±0.0001 cm- 3 , r c2 = 0.6547 ± 0.0004 arcmin, 
= 1.5474 ± 0.0001, xled = °- 71 > dof=7. And the best-fit 
profile is shown as the solid line in Fig. 4. 

4.2. Total mass 

Once we have derived the radial temperature profile T(r) 
and electron density profile n e (r) for RXCJ2228.6±2036, 
the integrated total mass of this cluster at radius r can 
be calculated under the assumptions of hydrostatic equi- 
librium and spherical symmetry by the following equation 
(Fabricant et al. 1980): 



M tot {< r) 



k B Tr 2 
G\xm v 



d{\nn e ) d(lnT) 



dr 



dr 



(9) 
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Table 3. The best-fit parameters for the central region of RXCJ2228.6+2036 by Model 2 . The lowT is fixed on 0.01 keV and 
1' = 332.7 kpc. The errors are at the 90% confidence level. 



T -^niefcol 


lowT c/ 


highT c/ 


A 


norm 


M 


An H 


xLd/dof 


(keV) 


(keV) 


(keV) 


(solar) 


(10~ 3 cm~ 5 ) 


(M ) 


(10 22 cm- 2 ) 




0' - 0.5' 8.54+^;^ 


0.01 (fix) 


= Tmekal 


0.32 ±0.19 


0.70t^ 


14 +ti4 ' u 


O.O(fix) 


1.12/250 




12 5 12 

channel energy (keV) channel energy (keV) 



Fig. 3. The spectra of the central region (r < 0.5') for joint fit of pn (bold crosses) and MOS (faint crosses) of RXCJ2228. 6+2036. 
a) fitted by a single-temperature model; b) fitted by a cooling flow model with an isothermal Mekal component. In b) we plot 
the isothermal and the cooling flow components respectively, and the lower lines below the crosses represent the multiphase 
components of pn (bold line) and MOS (faint lines), which show that the multiphase components only contribute a little to the 
emissions. 




ties of all the other properties of RXCJ2228. 6+2036 are 
also calculated from the 250 simulated clusters. 

Then, we derived the total mass profile of RXCJ2228.6 
+2036, shown in Fig. 5, and M tot = (1.36 ± 0.51) x 10 15 
M at 6' at the 68% confidence level. 



1 .00 r- 



0.10 - 



Radius (Mpc) 



Fig. 4. The deprojected electron density profile of 
RXCJ2228. 6+2036 after PSF correction. The error bars are 
at the 68% confidence level. The solid line is the best-fit profile 
from a double-/? model fitting. 



here ks is the Boltzmann constant, G is the gravitational 
constant and fi is the mean molecular weight in units of 
the proton mass m p (we assume fi = 0.6 in this work). 
The mass uncertainties are obtained from the uncertain- 
ties of the temperature and the electron density calculated 
by Monte-Carlo simulations. We obtained 250 redistribu- 
tions of the parameterized temperature and electron den- 
sity profiles by fitting to the data points varied within the 
Gaussian error bars of the measurements. The uncertain- 



0.01 r- 



0.1 



1.0 



Radius (Mpc) 



Fig. 5. The total mass profile of RXCJ2228.6+2036, and the 
error bars (dotted lines) are at the 68% confidence level. The 
vertical line indicates rsoo = 1.61 ± 0.16 Mpc. 

A physically meaningful radius for the mass mea- 
surement is defined as rsoo, the radius within which the 
mean gravitational mass density < p to t >— 500p c , where 
p c — 3H 2 /(8ttG) is the critical cosmic matter density. For 
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our calculations, we use the value of p c at the cluster red- 
shift, i.e., p c = 9.2 x 1CT 30 g cm~ 3 . This radius is still 
well covered by the observations. From the mass profile 
wc derive r 500 = 1.61 ± 0.16 Mpc for RXCJ2228.6+2036, 
corresponds to 4.8' ± 0.5', and the total mass within it 
is about M 500 = (1.19 ± 0.35) x 10 15 M Q . The mass 
derived from Chandra data by Morandi et al. (2007) is 
M tot = (4.90 ± 4.35) x 10 14 M Q within r 50 o = 1033 ± 464 
kpc, still consistent with ours within the error bars. In our 
analysis, the extrapolated value of r v » r — 2.48 ±0.38 Mpc 
= 7.5' ± 1.1' and M mr = (1.55 ± 0.72) x 10 15 M Q , which 
agrees with what Pointecouteau et al. (2002) derived. 

4.3. gas mass and gas mass fraction 

In galaxy clusters, gas is an important component involv- 
ing complex physics. From the electron density we calcu- 
late the gas mass and the gas mass fraction defined as 
fgas(r) = M gas (r)/M t ot(r). Fig. 6 shows the gas mass 
fraction profile of RXCJ2228.6+2036. The gas mass frac- 
tion at 6' is f gas = 0.17±0.06, consistent within the error 
bars with the results of LaRoque et al. (2006) from the 
Chandra X-ray and OVRO/BIMA interferometric SZ ef- 
fect measurements: f gas = 0.138 ± 0.009 from the X-ray 
data and f gas = 0.188^11 horn the SZ data. It also 
agrees with the WMAP measured baryon fraction of the 
Universe f b = 6 /O ro = 0.166 (Spergel et al. 2003). 



0.30 p 



0.20 



2 0.15 



0.05 



0.00 




0.5 1.0 1.5 

Radius (Mpc) 



2.0 



Fig. 6. The gas mass fraction profile of RXCJ2228. 6+2036. 
The dashed horizontal line indicates the WMAP measurement 
of fi, — Q. b /Q, m = 0.166 (Spergel et al. 2003) and the vertical 
line indicates r^oo = 1.61 ± 0.16 Mpc. The error bars (dotted 
lines) are at the 68% confidence level. 



5. Discussion 

5.1. PSF-corrected spectra 

The spatially resolved spectral analysis is affected by the 
PSF. To correct for the PSF effect, we first calculate the 



redistribution matrix, Fij, which is the fractional contri- 
bution in the ith ring coming from the jth ring (Pratt & 
Arnaud 2002). We get this redistribution from our best fit- 
ted double-/? model of electron density (converted to emis- 
sion measure profile) and the PSF matrices (Ghizzardi 
2001). Since we have divided our cluster into 5 regions, we 
can obtain the fractional contribution in each ring coming 
from all the bins, and in Fig. 7, we plot the contribution 
from the bin, all inner and outer bins. 



1 .0 



- 0.8 - 



c 0.6 - 



- 0.4 - 



0.0 



1 2 3 4 5 6 

Radius (arcmin) 

Fig. 7. Redistributions due to the XMM-Newton PSF: the di- 
amonds represent the contribution coming from the bin, the 
stars from the inner bins and the triangles from the outer bins. 



Here, d is the observed spectrum of the ith ring af- 
ter a double-background subtraction, Si is the spectrum 
without PSF effect, so for our cluster RXCJ2228.6+2036, 
we have: 
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(10) 



From this function we can derive Si, which indicates the 
spectra after the PSF correction. Then, we deproject these 
spectra Si by the deprojection technique described in Sect. 
2.3, and thus derive the PSF corrected deprojected spec- 
tra. Fitting these spectra with Models, we obtained the 
PSF-corrected deprojected temperatures, shown as stars 
in the upper panel of Fig. 2. 

We find that the PSF-corrected temperatures agree 
with the measurements without PSF-correction. This may 
be due to the broad width of the regions we chose. However 
it should be noted here that these spectral fits are not 
as good as those in Sect. 3.1 because the PSF-correction 
procedure introduces significant uncertainties, mainly due 
to the inversion process (see Eq. 11). 
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5.2. Gas pressure and comparison with the SZ data 

With the temperature profile T(r) and electron density 
profile n e (r), we can derive the gas pressure profile of 
RXCJ2228.6+2036 as: 



P(r) = n e (r)k B T(r). 



(11) 



In terms of observables, the gas pressure can be checked 
against the SZ effect (Sunyaev & Zel'dovich 1972) coming 
from the cluster. Indeed the SZ effect is directly propor- 
tional to the integrated pressure over the line of sight: 



Ssz(r) = / k B T(r)n e {r)dr ■ f sz {v,T) 



f P(r)dr ■ f sz (v,T), 



(12) 



where ks, rn e , c and ot are the Boltzmann constant, 
the electron mass, the speed of light and the Thomson 
cross-section; fsz{v,T) represents the SZ spectral shape 
(including the relativistic corrections as computed by 
Pointecouteau et al. 1998). 

Here we integrated the gas pressure of RXCJ2228.6 
+2036 using Eq. (13), convolved with the PSF of the 45 m 
radio telescope NRO (the beam size at 21 GHz: Ofwhm ~ 
80 arcsec) and then compared it with the SZ radial profile 
(Pointecouteau et al. 2002) in Fig. 8. The diamonds are 
from the SZ data and the solid line represents our result. 
We found a remarkable agreement within the error bars. 




7 .3 

Radius (arcmin) 



Fig. 8. The integrated X-ray pressure profile of RXCJ2228.6 
+2036 convolved with the PSF of the SZ telescope (the solid 
line) and compared with the SZ radial profile (the diamonds) 
derived by Pointecouteau et al. (2002). All the errors are at 
the 68% confidence level. 



The good agreement of the X-ray and SZ surface 
brightness profile in Fig. 8 may lend itself to test biases in 
the derivation of the pressure profile from the X-ray data. 
The most interesting aspects concern a bias in the temper- 
ature measurement in the presence of a multitemperature 
ICM (e.g. Mazzotta et al. 2004, Vikhlinin 2006) and the 



overestimate of the gas density due to the enhancement of 
the surface brightness if the gas is clumpy. In the following 
we will investigate how these two effects are modifying the 
comparison of the X-ray and SZ data. 

If we assume that the ICM is in rough pressure equi- 
librium the two bias effects on the temperature and the 
density are actually linked (for n e ■ T = constant). While 
local unresolved density inhomogeneities will lead to an 
overestimation of the density and the prediction of a too 
high SZ-signal, the temperature of a clumpy medium will 
be underestimated compared to a mass average and result 
in an underprediction of the SZ signal. So both effects are 
at least partly compensating each other in our study. 

Quantitatively the overestimation of the gas density is 
given by: 



c 



<nl> 



<n e > 



2 • 



(13) 



where the overestimation factor is C — \C' . To quan- 
tify the underestimation of the temperature for this hot 
cluster, we can approximately use the approach of 
Mazzotta et al. (2004) (Eq.14) which yields: 



T 



< n, 



X n 1 - 75 > 



< n 



2.75 



> 



(14) 



where T s i is a good approximation of the spec- 
troscopic temperature as would be obtained from 
data analysis of Chandra and XMM-Newton observa- 
tions for a multi-temperature plasma, and T the 
mass weighted mean. 

As an example we calculate these effects for a homo- 
geneous distribution with a lower and higher cutoff of, 
Ti =< T > -e and T 2 =< T > +e, respectively. A 
more general distribution can also be seen as a super- 
position of many of these 'top-hat' distributions. Fig. 9 
shows the enhancement factors C and R as a function 
of the distribution width parameter, e. We note that the 
two effects don't cancel each other, but the effect of the 
temperature underestimate is about 2-3 times larger than 
the overestimat-e due to dumpiness. Still the overall ef- 
fect is not dramatic and does therefore not provide a very 
good diagnostics. Even for a broad temperature distribu- 
tion with e/ < T >« 0.75 for example, covering a tem- 
perature range (from the lower temperature to the higher 
temperature) of a factor of 7, we obtain an SZE underes- 
timate of about 30% and a gas mass overestimate of 
about 18%. 

This has also implications on the mass measurement. 
While the pressure profile and its derivative can be directly 
taken from the SZ-profile, we still require an independent 
absolute temperature measurement for the normalization 
of the mass profile. The above calculation shows now, that 
we don't obtain very precise new information on a possibly 
low biased temperature due to a multiphase ICM from 
having simultaneous X-ray and SZ observations. In the 
above example a temperature and mass underestimation 
of 40% is only indicated by an SZ deviation of 30%. 
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Fig. 9. Overestimation factor of the gas density C and underes- 
timation of the spectroscopic like temperature versus the mass 
weighted temperature R as a function of the width of a homo- 
geneous temperature distribution in the presence of pressure 
equilibrium. The combined effect C x R is for the underesti- 
mation of SZE. For the definition of the parameters see the 
text. 



5.3. Gas entropy and the S — T relation 

Following Ponman et al. (1999), we denned the entropy of 
the gas in clusters as: 



S(r) 



n e (r)3 



(15) 



This entropy corresponds to the heat supplied per particle 
for a given reference density. Fig. 10 shows the entropy 
distribution as a function of radius, where the diamonds 
represent the entropy obtained from the spectra fitting 
results and the solid line from the best fitted T(r) and 
n e (r) profiles. 




Radius (Mpc) 

Fig. 10. The entropy distribution of RXCJ2228. 6+2036. The 
diamonds represent the entropy derived from the spectral fit- 
ting results and the solid line from the best fitted T(r) and 
n e (r) profiles. The vertical line shows rsoo = 1.61 ± 0.16 Mpc. 
The error bars are at the 68% confidence level. 



Pratt et al. (2006) have shown for the S — T rela- 
tion measured from a sample of 10 local and relaxed 
clusters observed by XMM-Newton, that 6*0.3 °c T x s4 , 
where So. 3 means the entropy at 0.3r 20 o and Tx is the 
mean temperature in the region of 0.1r 2 oo < r < 0.5r 20 o- 
For RXCJ2228.6+2036, So.3r200 = 959 ± 130 keV cm 2 



and T x = 8.91 



+ 1.91 
-1.33 



keV. The So.3r200 versus Tx for 



RXCJ2228.6+2036 is plotted on the S- T relation derived 
by Pratt et al. (2006), shown in Fig. 11. The diamonds and 
the best fit line (the solid line) are from Pratt et al. (2006), 
and the star indicates the result of RXCJ2228. 6+2036. 
It shows that our entropy value for RXCJ2228. 6+2036 
is consistent (within the la error bars) with Pratt et al. 
(2006) S — T relation at 0.3r 20 o, once corrected for the 
expected evolution in a self-similar scenario of structure 
formation. 




Tx (keV) 

Fig. 11. Comparison of the present result with the S — T re- 
lation of Pratt et al. (2006). The star indicates the result of 
RXCJ2228.6+2036, and the diamonds and the best-fitted S-T 
relation line (the solid line) come from Pratt et al. (2006). Here 
h(z) = [0.3(1 + zf + 0.7] 5. The error bars of the star are at 
the 68% confidence level. 



5.4. M - T and L - T relations 

From the above analysis, we derived the tempera- 
ture, mass and X-ray luminosity (see Table 1) of 
RXCJ2228.6+2036: within r 500 = 4.8', M 500 = (1.19 ± 
0.35) x 10 15 Mq, T500 = 8.92±1;£§ keV and L boli5Q0 = 
28.83±|;|| x 10 44 erg s" 1 . Therefore we can compare 
RXCJ2228. 6+2036 to the empirical scaling relations for 
massive galaxy clusters, e.g. M500-T500 and L bo i t500 -T 5Q0 
derived from XMM-Newton data, such as in Kotov & 
Vikhlinin (2005) based on 10 clusters at z = 0.4 - 0.7, 
Arnaud et al. (2005) based on 10 nearby clusters (z < 
0.15) and Zhang et al. (2008) based on 37 LoCuSS clusters 
at z ~ 0.2. The comparison is shown in Figs. 12-13. The 
diamonds and the solid line are from Kotov & Vikhlinin 
(2005), the triangles and the dashed line from Arnaud et 
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al. (2005), the dotted line from Zhang et al. (2008), and 
the star indicates the result of RXCJ2228.6+2036. 

It shows that our result is consistent with any of these 
previous studies within the scatter of the relations. The 
agreement of our L — T relation with that of Kotov & 
Vikhlinin (2005, with objects in the same redshift range as 
ours) is remarkable, in particular as X-ray luminosity with 
its square dependence on density is a parameter which is 
very sensitive to morphological disturbances and which 
thus generally shows a large scatter. 




10 

T (keV) 



Fig. 12. Comparison of the present result with literature M—T 
relations. The star indicates the result of RXCJ2228. 6+2036, 
the diamonds and their best-fitted M — T relation line (the 
solid line) come from Kotov & Vikhlinin (2005), the triangles 
and the dashed line from Arnaud et al. (2005), and the dotted 
line from Zhang et al. (2008). Here h(z) = [0.3(1 + zf + 0.7] 3 . 
The error bars of the star are at the 68% confidence level. 



5.5. M - Y relation 

The integrated SZ flux Ysz <x J ksTnedV oc M gas T, and 
thus the simplest X-ray analog is defined as Yx — M gas T. 
Kravtsov et al. (2006) show that Yx is the best mass proxy 
with a remarkably low scatter and the M — Yx relation is 
close to the self-similar prediction. 

For RXCJ2228.6+2036, M 500 = (1.19 ± 0.35) x 10 15 
M and Y x>mo = 21.84^;™ x 10 14 M keV. We plot 
M 5no versus Yx in Fig. 14 (shown as a star) and compare 
it with the M -Y x relations of Zhang et al. (2008) (the 
solid line), Kravtsov et al. (2006) (the dash-dotted line), 
Nagai et al. (2007) (the dashed line) and Arnaud et al. 
(2007) (the dotted line), which shows a good consistency. 

6. Conclusion 

We presented a detailed analysis of the XMM-Newton ob- 
servations of the distant galaxy cluster RXCJ2228. 6+2036 
(z = 0.421) using our deprojection technique. Through 
the spectral fitting we derived the deprojected tempera- 
ture profile T(r). Weighted by normalizations, we derived 
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1 I , , , , , i , , , , , , , , i , , , 
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L_bol * (1 +zY" (10" erg/s) 

Fig. 13. Comparison of the present result with the L — T rela- 
tion of Kotov & Vikhlinin (2005). The star indicates the result 
of RXCJ2228.6+2036, and the diamonds and the best-fitted 
L — T relation line (the solid line) come from Kotov & Vikhlinin 
(2005). The error bars of the star are at the 68% confidence 
level. 



, 10- 




10 

Yx « h(z) (10" Mo keV) 

Fig. 14. Comparison of the present result with the M — Y re- 
lation of Zhang et al. (2008) (the solid line), Kravtsov et al. 
(2006) (the dash-dotted line), Nagai et al. (2007) (the dashed 
line) and Arnaud et al. (2007) (the dotted line). The star in- 
dicates the result of RXCJ2228. 6+2036 with the errors of 68% 
confidence level. h(z) = [0.3(1 + zf + 0.7] i. 

a mean temperature within r 50 o, T500 = 8.92l};32 keV, 
which confirms within the error bars the previous results 
by Pointecouteau et al. (2002) and LaRoque et al. (2006). 

Then we calculated the cooling time of this cluster 
and obtained a cooling radius of 147+10 kpc. Fitted by a 
cooling flow model with an isothermal Mekal component, 
we derived the mass deposition rate ~ i2.OI.f2io ^eY ^ 1 
within r coo/ . 

Using the radial density profile n e (r) and radial tem- 
perature profile T(r), we obtained the mass distribution 
of RXCJ2228.6+2036. At r 500 = 1.61 ± 0.16 Mpc, the to- 
tal mass is M 500 = (1.19 ± 0.35) x 10 15 M , in agreement 
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with the results of Pointecouteau et al. (2002), derived 
from a combined SZ/X-ray spatial analysis, and the gas 
mass fraction is f gas = 0.165 ± 0.045. 

We discussed the PSF-correction effect on the spectral 
analysis and found that the PSF-corrected temperatures 
are consistent with those without PSF correction. 

We found a remarkable agreement within the error 
bars between our X-ray results and the SZ measurements 
in Pointecouteau et al. (2002), which is of prime impor- 
tance for the future SZ survey. The X-ray total mass and 
X-ray observables for RXCJ2228. 6+2036 closely obey the 
empirical scaling relations found in general massive galaxy 
clusters, e.g. the S-T, M-T, L-T and M-Y relations, af- 
ter accounting for self-similar evolution. 
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